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Results of a detailed examination of the cascade nucleosynthesis resulting from the putative 
hadronic decay, evaporation, or annihilation of a primordial relic during the Big Bang nucleosynthesis 
(BBN) era are presented. It is found that injection of energetic nucleons around cosmic time lO^sec 
may lead to an observationally favored reduction of the primordial ^Li/H yield by a factor 2 — 3. 
Moreover, such sources also generically predict the production of the ®Li isotope with magnitude 
close to the as yet unexplained high ^Li abundances in low-metallicity stars. The simplest of these 
models operate at fractional contribution to the baryon density Q,bh^ ^ 0.025, slightly larger than 
that inferred from standard BBN. Though further study is required, such sources, as for example 
due to the decay of the next-to-lightest supersymmetric particle into GeV gravitinos or the decay 



of an unstable gravitino in the TeV range of abundance fig/i 
both the ®Li and ^Li abundances in low metallicity stars. 



5 X 10 show promise to explain 



Big Bang nucleosynthesis has since long been known as 
one of the most precise and furthest back-reaching probes 
of cosmic conditions and the cosmic matter content of 
the early Universe. It thus, for example, has significantly 
contributed to the notion that a large fraction of matter 
in present-day galaxies is believed to be of non-baryonic 
nature as well as considerably limited some extensions 
of the standard model of particle physics. Paramount 
to having become such a useful tool of cosmology was 
and is, not only the examination of early synthesis of 
light elements in it's simplest standard version, but also 
in a variety of non-standard, alternative scenarios, relax- 
ing a priori non- verified assumptions entering the calcu- 
lations of a standard BBN (SBBN). Depending on the 
light element yields obtained in these latter scenarios, 
such calculations may then either favor a modified ver- 
sion of BBN, or strengthen the case for the SBBN. In ei- 
ther case, calculations of non-standard BBN may be used 
to place limits on the cosmic condition in the early Uni- 
verse. In this spirit, technically advanced calculations of 
BBN including decaying particles during, or after BBN, 
an inhomogeneous baryon distribution, small-scale anti- 
matter domains, neutrino degeneracy, or varying funda- 
mental constants (for reviews cf. Q), among others, have 
been performed over the years, rendering SBBN (also by 
the principle of Occam's razor) as our preferred scenario 
for BBN 

On the observational side significant advances have 
been made with the advent of high-resolution spectro- 
graphs on the Keck- and VLT-telescopes and the result- 
ing capability to perform D/H abundance determinations 
of unprecedented accuracy in a few simple high-redshift 
quasar absorption line (QAL) systems. Furthermore an 
independent determination of the fractional contribution 
of baryons to the critical density, fif,, from precision mea- 
surements of the cosmic microwave background radiation 
(CMB) by various balloon missions and the WMAP 2] 
sattelite has become feasible. These observations to- 
gether with the ever-continuing observational and the- 
oretical efforts to deduce precise primordial ''He/H- and 
^Li/H- ratios, leave BBN in an essentially observationally 



overconstrained state, opening the possibility to question 
internal consistency of the predictions of SBBN. If such a 
check is performed it shows that SBBN is to first approx- 
imation internally consistent, though inconsistencies or 
tensions between predicted- and observationally inferred- 
abundances may exist at higher order. It is the subject 
of this paper to propose a scenario of non-standard BBN 
which may remove some of these tensions. 



If the central value of the observationally determined 
D/H= 2.78toit X 10"^ by the average of five QAL 

systems is taken, SBBN predicts "^Li/Hw 4.16 x 10"^°, 
a "^He mass fraction Yp « 0.2480 at an f^b/i^ w 0.0218. 
I refrain from the common practice to give error bars 
due to nuclear reaction rate uncertainties on the theo- 
retical predictions, as there exist surprisingly large dif- 
ferences between the central values obtained by different 
groups, often larger than the quoted error bars, partic- 
ularly in the case of ^Li [|. The predicted ^Li should 
be compared to the observationally inferred primordial 
^Li/H= 1.23l^ jg X 10-1° [13 abundance from the in- 
ferred ^Li in atmospheres of extreme Pop II stars be- 
longing to the Spite plateau (where the analysis corrects 
for ^Li production by cosmic rays), or to the inferred 
^Li = 2.19lE5 j? ^ P] from low-metaUicity stars 

within the globular cluster NGC 6397, indicating if taken 
at face value, that there may be a problem in SBBN. 
This discrepancy may not be resolved by nuclear reaction 
rate uncertainties in the main lithium-producing reaction 
3He(a,7)^Be and only very unlikely due to uncertain- 
ties in the lithium-destroying reaction '^Be(d,p)2'*He 8]. 
It is conceivable, however, that it is due to other sys- 
tematic uncertainties entering the inference of primordial 
^Li abundances, such as stellar astration of ^Li in low- 
metallicity stars, or imprecise determinations of stellar 
surface temperatures in these stars. Recent 3D non-LTE 
(i.e. local-thermodynamic-equilibrium) calculations 
of lithium lines in low-metallicity stars indicate, however, 
that simplifying assumptions concerning the stellar atmo- 
sphere probably do not introduce excessively large sys- 
tematic errors. 
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Stellar '^Li may be destroyed by ^Li(p, a)^He when 
transported deep enough into the interior of the star, 
as for example due to stellar rotation induced turbu- 
lence or diffusion. The mixing possibilities proposed 
are numerous, however, many fail on the requirement 
to deplete ^Li in different stars of different surface tem- 
peratures and rotation velocities by a substantial factor 
(~ 2 — 4 = 0.3 — 0.6 dex) neither introducing scatter in 
the observationally inferred ^Li/H nor a trend of deple- 
tion with surface temperature. Moreover, ''Li has to be 
depleted without substantial depletion of the more frag- 
ile ^Li 0, 01 isotope observed in at least two of these 
stars along with that of ^Li. The large abundance of the 
^Li isotope at such low metallicity is anyway theoreti- 
cally challenging, even in the absence of stellar depletion 
(see below). Recent detailed stellar studies (15iil6j and 
an analysis ^3] of the samples of ^Li/H abundances em- 
ployed to infer the primordial ^Li/H ratio claim that, 
under certain conditions, a uniform depletion of 0.2-0.3 
dex may be possible, though there remains lack of obser- 
vational indication for this, in fact, to be the case 

Concerning the abundance of ^He there exists as well 
a mismatch between the relatively high prediction and 
the lower observationally inferred value Yp — 0.02390 ± 
0.0020 [11 (with another group finding Yp = 0.2421 ± 
0.0021 [iyl)- I will not much further elaborate on this 
discrepancy here, since the abundance of ''He, whose ob- 
servational determination is plagued by various system- 
atic uncertainties, will hardly change in the scenarios I 
consider. Last but not least, there is the vital comparison 
between the fib/i^as inferred from the CMB and fifc/i^as 
preferred by SBBN and D/H. The inferred baryon density 
when the WMAP data only is taken 2J (with a power-law 
ACDM model) is ^bh'^^ 0.024 ± 0.001, whereas a run- 
ning spectral index ACDM model results in a value lower 
by 0.001. When data from other experiments, of CMB 
and/or large-scale structure, probing smaller scales than 
WMAP are included, the estimated rif,/i^drops to around 
0.0225 (for combined likelihoods between fib/i^estimates 
from CMB and SBBN cf. to Ref. (3). In general the 
agreement between CMB and SBBN (plus D/H) is ex- 
cellent, with nevertheless CMB inferred r^b/i^values typi- 
cally staying on the high side of those preferred by SBBN. 
Thus non-standard BBN scenarios are not allowed to op- 
erate with ilf,/i^much different than those of SBBN. 

At such large flbh^ihe bulk of ^Li is produced as ''Be 
(which is after BBN converted to "^Li). It is a not widely 
recognized effect that ''Be may be prematurely destroyed 
via the reaction chain ^Be(n,p)^Li and ^Li(p, a)^He by 
significant factors when towards the end of BBN, at ap- 
proximate temperatures 30 — 50 keV a small excess of free 
neutrons n/p ~ 10~^ as compared to SBBN exists. This 
phenomena has been observed in both studies of inho- 
mogeneous BBN [2p| as well as studies with hadronically 
decaying particles [2l|. Such an excess could thus pos- 
sibly result from the hadronic decay or evaporation of a 
relic of the early Universe. Thermal neutrons injected at 
these temperatures are predominantly incorporated into 




FIG. 1: Light element evolution as function of temperature 
with (solid) and without (dashed) a thermal neutron source. 
Shown are, from top to bottom (at lower temperatures), the 
number ratios nn/np, nan^/rip, rin/np, m^^/up, and m-^^/np, 
respectively. See text for further detail. 



D via p(n,7)D, with a small fraction ^ 10^^ — 10^^ 
causing the premature conversion of ''Be to ''Li. It is 
straightforward to show that to convert one ^Be nucleus 
{o'v)pn/{<7v)7Ben X H/''Be ueutrous are required. Due 
to the ^Be(n,p)^Li reaction being between mirror nuclei, 
the rate ratio in the above is very small w 1.4 x 10^^. 
This implies that the conversion of essentially all ''Be to 
^Li is accompanied by only a mild (^Be/H -independent) 
excess of D/H « 1.4 x 10^^, when compared to SBBN. Of 
course, not all of the thus produced ^Li will subsequently 
be destroyed via ''Li(p, a)*He, in particular at low tem- 
peratures where the Coulomb barrier is preventing the 
reaction. Neutrons which are injected at slightly higher 
temperatures are further processed into ^He. Neverthe- 
less, this possible excess in ''He is essentially negligible 
due to the small numbers of neutrons required to destroy 
^Be. 

I have slightly modified the Kawano code in order to 
test for this effect. In particular, I have injected ther- 
mal neutrons with rate dun/dt ~ exp(— t/T)/T and decay 
time T = 700 s employing a total injected neutron abun- 
dance corresponding to rj„/i^ « 10~^. In Fig. 1 the light- 
element synthesis in the presence (and absence) of such 
a neutron source is shown, assuming fl^h'^ = 0.026. It is 
seen that the neutrons in fact yield the desired destruc- 
tion of ^Bc and sonic comparatively smaller enhancement 
in ''Li. The yields for D/H and ^Li/H, with (and with- 
out) extra neutrons are 3.25 x 10~^, (2.09 x 10~^), and 
1.73 X 10^1", (5.92 X 10~i°), respectively, resulting in a 
0.53 dex "depletion" of the ^Li yield. The ''He abundance 
remains virtually unchanged (except for an increase by 
AYp « 0.002 on account of the increased f2b/i^with re- 
spect to the above quoted value). 

The effect is encouraging and may present a possible 
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resolution to the ^Li problem in SBBN. Nevertheless, 
many viable hadronically decaying, evaporating, or anni- 
hilating candidates inject energetic hadrons rather than 
thermal neutrons. A massive ~ 200 GeV decaying 
particle, for example, would lead to the injection of pro- 
tons and neutrons with typical energies ~ 8 GeV reaching 
up to several tens of GeVs. These are accompanied, of 
course, by much larger numbers of pions, e^, neutrinos, 
and photons. Cascade nucleosynthesis due to energetic 
electromagnetically interacting particles may not be of 
immediate interest here, not only because it is opera- 
tive at lower T ^ 3 keV, but also since the small required 
neutron densities to resolve the ^Li discrepancy may not 
imply much of an effect. Similarly, n^p interconver- 
sion by mesons |^ |22j , potentially important at higher 
temperatures, has a negligible effect due to the small- 
ness of the assumed perturbation. I caution though that 
these conclusions are dependent on the injected meson- 
and baryon- multiplicities, and the ratio of electromag- 
netically to hadronically interacting particles, and are 
strictly only valid in the context of strongly interacting 
jet dynamics. In contrast, spallation of ^He by energetic 
nucleons is important, and has heretofore only been con- 
sidered after BBN by the pioneering study of Ref. 
(cf. also to [131 ), and with an amount of injected nu- 
cleons far larger than of interest here. Similar holds for 
nonthermal nucleosynthesis, in particular, the reactions 
'^H(a,n)^Li and ^He(a,p)^Li with energy threshold in- 
duced by energetic '^H and ^He (themselves produced by 
the spallation of "^He) are of paramount importance for 
the ^Li abundance. 

The large abundance of ^Li in low-metallicity Pop II 
stars is a bit of a mystery. It has been observed in 
at least three stars at low metallicity [Fe/H]< —2 with 
the best case probably the star HD84937 with ^LifU^ 
0.052 ±0.019 [23, as well as in two stars at higher metal- 
licity [Fe/H]ss —0.6 jl^, with all abundance ratios co- 
incidentally being similar '27^. ^Li is traditionally not 
thought to be of primordial origin due to the cross sec- 
tion D(Q!,7)^Li, absent of threshold, being small [28j . 
^Li is thus believed to have its origin in cosmic-ray nu- 
cleosynthesis, being produced along-side with ^Li (to be 
added to the BBN yield), ^Be, ^°B, and "B via spallation 
(p, a-l-CNO LiBeB) and fusion {a + a ^ Li) reactions 
by galactic cosmic rays. It is rather controversial as to 
whether , or not [sO, E| , it is possible to produce the 
large observed ^Li abundance in Pop II stars via cosmic 
rays generated by thermonuclear or core-collapse super- 
novae, even when so far unknown cosmic ray populations 
(e.g., low-energy and metallicity-dependent) are postu- 
lated. The problem, in general, seems to be that in order 
to reproduce the approximate linear relationship between 
^Be/H and [Fe] metal-enriched cosmic rays are strongly 
favored. Such metal-enriched cosmic rays yield typical 
spallation ^Li/^Be-ratios ~ 5 — 20 agreeing with the ra- 
tio w 6 in the solar system, but not with the ratio ~ 80 
as observed in Pop II stars. It may be that there exists 
an enhanced fusion contribution (with cosmic rays at ~ 



30 MeV/nucleon) at high redshifts, and though it may 
be energetically problematic when associated with stel- 
lar evolution [30|, it could be due to viriahzation shocks 
during structure formation Nevertheless, the ori- 

gin of ®Li in Pop II stars seems currently controversial. 
This situation has led me to consider alternative pro- 
duction of ^Li due to non-thermal nucleosynthesis after 
BBN 33, 34|, resulting from electromagnetic cascade nu- 
cleosynthesis. 

Production of ^Li due to nuclear spallation and fu- 
sion reactions towards the end of BBN may as well 
be very efficient. I have performed a detailed Monte- 
Carlo analysis of the cascade nucleosynthesis resulting 
from the decay, evaporation, or annihilation of strongly 
interacting particles or defects. To reach a precision 
to allow for a meaningful comparison to the compar- 
atively accurate observational determinations of D/H, 
^Li/H, ^Li/H, and ilbh^i care had to be taken. The 
analysis includes a detailed modeling of elastic nucleon- 
nucleon and nucleon-"'He scattering. It includes also 
their important inelastic {pp,np — > pp,pn,nn -\- n's and 
p, n^He ■^He, "^H, D s -I- n's,p's + tt's) counterparts em- 
ploying detailed cross section data. A careful treatment 
of the recoil energies of ''lie (from elastic scattering), as 
well as and '^He (from "'He spallation) has been incor- 
porated. These recoil energy distribution functions are 
important for a determination of ^Li yields resulting from 
the reactions ^R{a,nfU, ^lle{a,p)^U, '^lle{a,2pfRe, 
'*IIe(a,pn)^Li. Finally, as thermalization of nucleons, 
and nuclei, is a competition between nuclear scattering 
and Coulomb interactions with the plasma, in the case of 
protons and nuclei, and nuclear scattering and magnetic 
moment scattering on e^, in the case of neutrons, these 
latter interactions also had to be included [s^. A de- 
tailed account of the analysis, which is beyond the scope 
of the present paper, will be presented elsewhere. These 
processes have been cougled to the Kawano code [s^ (up- 
dated by the NACRE |37| reaction compilation) . Finally, 
the primary spectrum of injected nucleons is computed 
with the help of the jet fragmentation code PYTHIA [s^ 

The development of a nuclear cascade may be sum- 
marized as follows. For temperatures above T ^ 20 keV 
protons are essentially exclusively thermalized by electro- 
magnetic interactions with the still abundant pairs. 
At lower temperatures, thermalization for energetic pro- 
tons (and nuclei) iJkin ^ 1 GeV occurs via scattering on 
plasma protons and ''He, as well as Thomson scatter- 
ing on CMBR photons, while less energetic protons loose 
energy via Coulomb scattering and plasmon excitation. 
In contrast, neutrons thermalize on plasma protons and 
"'He up to temperatures of T « 50 — 60 keV. For higher T 
the dominant energy loss is due to ne^-scattering, with 
the most rapid fractional energy loss for energetic neu- 
trons (7 3> 1). Thermalization of nucleons and nuclei 
occurs rapidly when compared to the cosmic expansion 
rate. Interactions of rapid nucleons on protons (and ^He) 
at i?kin ^ 1 GeV are to a large fraction ('^ 0.75) inelastic, 
accompanied by the production of pions. During such 
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FIG. 2: Abundance yields of D/H, ^Li/H, and ''Li/'^Li in 
an 0,i,h^ — 0.026 Universe as function of the hadronic de- 
cay time r of a putative primordial relic. The models are 
decay of a m-)^ = 10 GeV particle (long-dashed), decay of a 
rrix = 200 GeV particle (solid), decay of a = 4TeV par- 
ticle (dashed-dotted), injection of monoenergetic nucleons of 
Ekin = 250 MeV (short-dashed), and extended power- law in- 
jection due to a = 200 GeV particle (dotted). Also shown 
are the two-sigma ran ges of the inferred primordial D /H and 
^Li/H abundances |3. llOl as well as the ^Li/^Li ratio as in- 
ferred in the low-metallicity star HD84937 12511 . See text for 
further details. 



scatterings an interconversion of protons to neutrons oc- 
curs frequently, such that energetic protons produce sec- 
ondary neutrons. For example, though the decay of a 
200 GeV particle generates only about « 1 neutron per 
annihilation, around « 1,0.6 secondary neutrons result 
at T « 20,40keV, respectively J39j, and « 3.5 asymptot- 
ically at low temperatures T ^ 0.1 — 1 keV. Here at higher 
temperatures the number of secondary neutrons reduces 
due to the rapid Coulomb losses of protons. Neutrons, 
on the other hand, do not possess a significant bias to- 
wards producing secondary neutrons in np inelastic inter- 



actions. Excess neutrons at T w 40keV are mostly due 
to inelastic processes on ^He, accompanied by the pro- 
duction of D and'^He (i.e. n-l-^He D+p+2n, ...), with 
a comparatively smaller amount of neutrons removed in 
pionic fusion processes (i.e. np Dtt*^, ...). One thus 
obtains approximately a ratio n/Dw 3.6 for a 200 GeV 
particle at T « 40keV, with similar ratios for n/^H and 
n/^He. As the '^H and '^He are energetic they may yield 
the production of ^Li. Nevertheless, ^Li production (and 
survival) may only be efficient at somewhat lower temper- 
atures. Due to Coulomb losses of energetic and '^He 
production is only efficient at T ^ 20 keV, whereas sur- 
vival of the freshly synthesized ^Li against destruction via 
^Li(p, a)'^He is only nearly complete for r< lOkeV. The 
production of ^Li at temperatures T « 10 — 20 keV for a 
200 GeV particle is found to be approximately 2 x 10"'' 
per decaying particle, becoming significantly lower at 
lower temperatures (e.g. 3 x 10~^ at T w 1 keV). Cascade 
yields are subject to some nuclear physics data uncertain- 
ties which in the case of ^Li may be of the order of a factor 
two. In particular, it may be that ^Li yields are under- 
estimated due to an experimentally incomplete determi- 
nation of the high-energy tail of the energy distribution 
of energetic '^H and '^He produced in ^He spallation. 

The developed code allows me to present detailed pre- 
dictions on the BBN in the presence of decaying parti- 
cles. Figure 2 shows the light-element yields for a variety 
of decaying particles as a function of particle life time 
T. The panels show, from top-to-bottom, final abun- 
dances of D/H, ''Li/H, and ^Li/''Li, with the understand- 
ing that Yp is virtually unchanged when compared to 
SBBN at the same flbh'^. In all models f^b/i^^ 0.026 
has been assumed. Hadronically decaying particle yields 
(with the simplifying assumption that x ~* 99 yields the 
production of a pair of quarks, the up-quark for definit- 
ness) are shown for three particle masses: = 10 GeV 
with n^h"^ = 7.5 X 10~^ (long-dashed), = 200 GeV 
with n^h'^ = 1 X 10""* (solid), and ^ 4TeV 
with il^/i^ 6 X 10^"* (dashed-dottcd). It is evident 
that for decay times around r « lO'^s an efficient de- 
struction of ^Li is obtained. For r much shorter than 
lO'^s the destroyed ^Be is regenerated, whereas for r 
much longer, incomplete ^Li burning in the reaction chain 
^Be(n,p)^Li(p, Q;)^He results in only partial reduction of 
the total ''Li yield. As anticipated, the destruction of ^Li 
is accompanied by production of D. When compared to 
the injection of thermal neutrons, D/H yields are higher. 
This is due to D generated in the nuclear cascade it- 
self (i.e. by ''He spallation and pionic fusion). Cascade 
generated deuterium (as well as '^H, •^He, and ®Li) is sub- 
stantially reduced per injected neutron for sources which 
inject nucleons with a soft spectrum. For example, I have 
also employed a soft source with monoenergetic nucleons 
of 250 MeV. Results for this case are shown by the short- 
dashed line, assuming n^h'^/m^ « 7.5 X lO-^'GeV-i and 
the injection of one np pair per decay ^ cascade 

n/D « 10 ratio at T « 40keV is obtained in such scenar- 
ios. The more pronounced depth of the ^Li dip in Fig. 
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2, by comparable (or smaller) D production, indicates 
that soft sources may affect a larger destruction factor in 
^Li while respecting an appropriate upper limit on the 
primordial D/H-ratio. 

Concerning production of ®Li it is observed that yields 
of the order of ^Li/^Li « 5 x 10^^, as observed in extreme 
low-metallicity Pop II stars, are obtained for decay times 
in the range t « 1 — 2 x lO'^s. For such t, substan- 
tial ^Li destruction is also observed, since minimal ^Li 
production usually occurs for t between 800 and lO'^s. 
The slight mismatch between the "perfect" decay time 
for ^Li destruction when compared to that for ^Li pro- 
duction is related to the more efhcient destruction of ^Li 
(via ®Li(p, a)^lle) as compared to the analogous reaction 
for ^Li. It also depends on the total ^Li yield per decay- 
ing particle and may disappear if their exists a factor two 
underestimate in this quantity. However, it is not neces- 
sary that both decay times are equal, as even when not, 
observationally acceptable scenarios may result. This is 
illustrated in Table 1, which shows the abundance yields 
of a few selected models, generally in good agreement 
with the observational value of all three isotopes D, ^Li, 
and ^Li. If it all than the abundance of D ~ 4 x 10^^ in 
these decaying particle scenarios is somewhat large. 

It is tempting to adjust for this slight mismatch by re- 
sorting to somewhat less in the literature discussed possi- 
bilities of hadronic energy injection than that of a massive 
decaying particle (such as the gravitino). These may in- 
clude unstable Q-balls :42j or semistable strange-quark 
matter nuggets formed during a QCD-transition [43l |. 
among others. What is desired is a prolongation of the 
injection, as compared to the decay of a particle. This 
may be, for example, due to a population of Q-balls 
of varying sizes. As a detailed analysis is beyond the 
scope of this exploratory paper, I simply model the tem- 
poral injection of such a putative source by dnq/dt — 
Uq /t mm[{t / t)^^ , 1], including a long time injection tail. 
Here the spectrum of the primary and secondary cascade 
products is taken to be that of a 200 GeV qq event and 
fig/i^ = 10~^ has been assumed. Results for this case 
are given by the dotted line in Fig. 2, and for one par- 
ticular T on the last row of Table 1. It is seen, that 
for such "extended" emission all, the D/H, ''Li/H, and 
^Li/^Li abundance constraints may be well satisfied. As 
another possibility I mention sources which have an un- 
usually large ratio between electromagnetic and hadronic 
energy injection (~ 10^ — 10^). If emission lasts to around 
times of lO^s, a substantial fraction of D (as well as ^Li) 
may be photodisintegrated. Nevertheless, such sources 
have subsequently to stop radiating fairly abruptly, as 
otherwise an overproduction of ^Li (either via ''Li photo- 
disintegration or electromagnetic cascade nucleosynthe- 
sis [S^) results. Note, that this is also the reason that 
the recently proposed ^Li photodisintegration as a solu- 
tion to the comparatively high predicted ^Li abundance 
may not apply ^j. 

Extended emission during BBN may result as well 
from annihilating particles. Since a supersymmetric neu- 



TABLE I: Light-element abundances for a few selected de- 
caying particle and extended emission scenarios. All models 
assumme Qbh^ = 0.026. See text for further details. 



model 


r 


(s) 


D/H 


^Li/H 


^Li/^Li 


thermal neutrons 


7- 


10^ 


3.3 ■ 10" 


5 


1.7- 10" 


10 


<io-3 


soft source 


1.8 


• 10^ 


4.0 • 10" 


5 


1.9 ■ 10" 


10 


0.03 


decay — 10 GeV 


1.5 


• 10^ 


5.1 • 10" 


5 


1.4 ■ 10" 


10 


0.04 


decay = 200 GeV 


1.5 


• 10^ 


4.1 ■ 10" 


5 


2.1 ■ 10" 


10 


0.03 


decay — 4 TeV 


1.5 


• 10^ 


4.0 • 10" 


5 


2.4 • 10" 


10 


0.03 


extended emission 


7.1 


• 10^ 


3.6 • 10" 


5 


1.5 • 10" 


10 


0.05 



tralino, in case it is the lightest supersymmetric particle 
(LSP) and when R-parity is conserved, leads to residual 
annihilation during BBN it is of interest to explore this 
case. In Fig. 3 I show the nucleosynthetic signatures of 
annihilating neutralinos in D/H, ""Li/H, and ^Li/^Li as a 
function of Vlbh? ■ Here the light solid, long-dashed, and 
dotted lines correspond to residual hadronic annihilation 
from light 5 and 10 GeV neutralinos with varying s-wave 
annihilation cross sections {av) = 3 and 5 x 10"^^cm^/s. 
These values of {av) correspond approximately to those 
required for a thermal freeze-out of neutralino annihi- 
lations before (and after) the QCD transition to yield 
Vl^h? ss 0.1126. Though such models formally violate the 
LEP lower bound on the mass of the lightest neutralino 
^ 50 GeV, this latter bound implicitly assumes gaug- 
ino mass unification at the GUT scale. When this as- 
sumption is dropped, my_ « 5 GeV neutralinos remain a 
viable dark matter candidate 45] . It is seen, that with re- 
spect to SBBN, the ^Li overproduction problem is hardly 
alleviated, while such models produce ^Li in excess of 
that observed in Pop II stars. ^Li overproduction results 
particularly for lower and larger {av), as here resid- 
ual annihilation is at it's strongest. To obtain the desired 
destruction of ^Li models would lead to a factor 10 — 20 
overproduction of ^Li. Models with p-wave annihilation 
{av) ~ Tjm^ fare somewhat better, but at the expense 
of introducing exessively large annihilation rates. In any 
case, residual annihilation of dark matter particles seems 
not to be capable of explaining at the same time, the 
apparent SBBN ^Li overproduction problem and the ^Li 
abundance in low-metallicity stars. 

I have so far implicitly assumed that ratios of D/H 
^ 3 — 4 X 10"*^ are observationally allowed. This may 
seem somewhat at odds with the observational determi- 
nation of D/H= 2.78;5J:38 X 10-5 H]. In this context, 
it should be noted (as Ref. did), that the formal er- 
rors given above may represent an underestimate of the 
true error (due to further unknown systematics). As a 
spurious trend of D /H with total hydrogen column den- 
sity of the Lyman-a absorbers exists, with lower col- 
umn density systems giving higher inferred D/H (such 
as Q1009-I-2956 with 4.0 ± 0.65 x lO'^ li^ or PKS 1937- 
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FIG. 3: Abundance yields of D/H, '^Li/H, and ''Li/'^Li 
in the presence of residual neutralino annihilation for a 
neutralino fractional contribution to the critical density of 
Qxh'^ = 0.1126. Three models are shown: = 5 GeV 

and (av) = 3 x 10~^®cm^/s (light solid), = 5 GeV and 
(av) = 5 X 10"^'^cm^/s (dashed), and = 10 GeV and 
(av) = 3 X 10"^''cm^/s (dotted). For comparison, the SBBN 
yields are shown by the heavy solid lines. Note the absence 
of appreciable ®Li production in SBBN. Observational data 
on the D/H, ^Li/¥L, and ^Li/^Li ratios indicated in the figure 
are as in Fig. 2. 



1009 with 3.25 ± 0.3 x IQ-^ [13), I regard primordial 
values of D/H 3 — 4 x 10~^ as viable. Similarly, a 
baryonic density parameter of fli,h^= 0.026 may seem 
compatible with the determination of ^bh^when WMAP 
data is only considered, but at odds with the estimate 
of the combined data of WMAP and large-scale struc- 
ture surveys (as well as Lyman-a forest surveys) which 
yields 0.0226 ± 0.0008. Nevertheless, since such analy- 
sis invokes assumptions of the underlying cosmological 
model, formal error bars should be taken with a grain 
of salt. Evidently, scenarios which are accompanied by 



additional late-time photodisintegration of D (and ^Li) 
may operate at lower fli,h^. 

It is interesting to speculate on which relic, in fact, 
could cause the depletion of ^Li and concommitant pro- 
duction of ^Li. Though residual annihilation of a dark 
matter particle usually predicts too large of a ^Li abun- 
dance, when requiring a factor 2 — 3 ^Li depletion, it may 
still represent the main source of ^Li in low-metallicity 
stars. Alternatively, an initially high ^Li may (along 
with the ^Li ) be somewhat depleted by stellar processes 
to reach agreement with observations on both isotopes. 
Other possibilities include the evaporation of Q-balls or 
strange quark matter nuggets, but will not be further 
discussed here. In light of the widely believed existence 
of a supersymmetric (SUSY) extension of the standard 
model of particle physics, particle decay may seem par- 
ticularly promising, as there exist a variety of possibil- 
ities involving the gravitino. When spontaneous SUSY 
breaking in a hidden sector is communicated to the vis- 
ible sector via gravitational interactions, or when break- 
ing occurs due to the super- Weyl anomaly, the gravitino 
is heavy ~ 0.1 — 100 TeV and long-lived. Moreover, it 
is easily produced in the primordial plasma due to two- 
body scatterings at high temperture Tr with abundance 
~ 2 x 10"'^(Tfl/10«GeV)(TeV/TOg.) 49j. It's de- 
cay into a gauge boson B and it's superpartner B oc- 
curs with life time r, 



G 



4 X 10^sN-^{m^/lTeV) 



where Nc is the number of available decay channels {Nc 
= 1 for the U{1)y gauge boson and Nc = 8 for gluons). 
For decay into gluons one finds r^, w 700 s for a 4 TeV 
gravitino |48j |. seeming somewhat large for the simplest 
models of supergravity and somewhat small for models 
of anomaly-mediated SUSY breaking. Other scenarios 
may result when the gravitino itself is the LSP (light- 
est supersymmetric particle), such as in gauge- mediated 
SUSY breaking (GMSB). In this case, the particle next 
in the mass hierarchy (NLSP) undergoes thermal freeze- 
out of annihilation reactions, with typical abundances 
of the order i^^h^ 10^'^ — 10"* independent of the 
reheat temperature. The NLSP later on decays into 
the gravitino LSP and it's superpartner with decay time 
« 240s(m(5/GeV)2(mx/300GeV)-5, thus making de- 
cay occur naturally at the desired lO^s for GeV graviti- 
nos. In GMSB scenarios the NLSP is frequently the bino 
B, right-handed stau f, and in a smaller part of the pa- 
rameter space the higgsino H. Typical annihilation rates 
for B, f, and H are in the ballpark 10^^^ for the for- 
mer and lO^^^cm'^/s for the latter two particles, yield- 
ing thermal freeze-out fJ^ of ~ 1 and 10~^, respectively. 
Nevertheless, a significant spread around these reference 
values exists. As I have shown an Bh^^h? « 2 x lO^'' 
is preferable, where Bh is the hadronic branching ratio 
capable of producing nucleons (i.e. excluding hadronic 
decay only generating pions as in the case of the r) . In 
the case of the f and the H the typical abundance comes 
close to this value. The f , however, has only a small 
phase-space supressed Bh ^ 5|/(327r^) ~ 10"'^ due to 
the coupling of r to Vr and , with the latter decay- 
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ing into qq in seventy per cent of all cases. Here (72 is the 
SU{2) gauge coupling constant. An NLSP stau may thus 
only work for Clf ^ 0.1 — 1, demanding it to be rather 
heavy m-f- ~ 1 TeV. For a Higgsino one has Bh ~ 0.1 — 1 
due to it's decay into heavy quarks or massive gauge 
bosons, such that a Higgsino NLSP may have the de- 
sired properties. Such scenarios may also be consistent 
with the gravitino being the dark matter particle and 
successful leptogenesis occurring at reheat temperatures 
Tr ~ IQiOGeV 153. Finally, a typical bino abundance 
seems too large since the bino is usually appreciable 
(~ 0.1) due to decay into the Z boson. Nevertheless, in 
this case Bh may be phase space suppressed for not 
much larger than the Z mass. Particularly interesting in 
this case of Bh ^ 1 is also the possible later photodis- 
integration of D, possibly removing the requirement to 
operate at ^l^h^ ^ 0.025. Such decay is also associated 
with an, albeit small, component of gravitino hot dark 
matter jsj. Last but not least, the particle with the re- 
quired properties may also be a super- WIMP occurring 
in Kaluza-Klein theories [s^. 

In summary, I have presented first results of a newly 
developed Monte Carlo code examining the nuclear cas- 
cade development and cascade nucleosynthesis resulting 
from the putative injection of energetic nucleons during 
the epoch of BBN. Such an injection may result by a va- 
riety of means, such as, the decay of gravitinos, the de- 
cay of NLSP's to gravitinos, the evaporation of Q-balls 
or strange quark matter nuggets, or the residual anni- 
hilation of light neutralinos, among others. This most 
detailed numerical tool of this sort to date, has been ap- 
plied in the regime of a " comparatively" weak perturba- 
tion during BBN leaving the primordial *He abundance 
virtually unchanged when compared to SBBN. However, 
sources of this kind may still have dramatic effects on the 
synthesized D, ^Li, and ^Li yields. In particular, I have 
found that due to injection of neutrons near r « lO'^s an 
efficient reduction of the final ^Li yield may result. Such 
a factor 2-3 reduction of ^Li is just what is needed to re- 



solve the tension between the observationally inferred low 
^Li/H ratios in low-metallicity Pop II stars and the theo- 
retically predicted high ^Li/H ratio in SBBN. Of course, 
this holds true only if the ^Li in these stars is not subject 
to some significant stellar depletion. In the simplest sce- 
narios, destruction of ^Li is found to be accompanied by 
some production of D, though in more complicated sce- 
narios such additional D may later be partially photodis- 
integrated. If photodisintegration is absent, viable sce- 
narios may only result for ^l^h^ ^ 0.025, such that future 
(and current) high precision determinations of flbh^hy 
measurements of CMB anisotropics are important for the 
evaluation of these models. Injection of energetic nucle- 
ons towards the end of BBN may also result in the syn- 
thesis of the ^Li isotope, with resultant ^Li/^Li isotope 
ratios comparable to those observed in low-metallicity 
Pop II stars. This may present a viable alternative to 
the problematic explanation of ^Li abundances in Pop II 
stars due to traditional cosmic ray nucleosynthesis.. It 
is intriguing to note that BBN with a weak non-thermal 
hadronic source, shows good potential to resolve two cur- 
rent discrepancies in nuclear astrophysics. In this con- 
text, particularly promising seems the decay of NLSP's 
to LSP gravitinos of mass ~ IGeV with NLSP freeze-out 
abundance either small « 2 x 10^'* when an appre- 
ciable hadronic (nucleonic) branching ratio Bh exists, or 
larger when Bh ^1. I have shown that this may the 
case for the NLSP being either the bino, stau, or Hig- 
gsino. Another viable possibility is the decay of ^ TeV 
gravitinos of abundance fl^ sa 5 x 10^* produced during 
reheating at a reheat temperature of Tr = lO^GeV. I be- 
lieve that such processes deserve further consideration as, 
if indeed they occurred, they may provide an invaluable 
source of information about the evolution of the early 
Universe and the properties of a primordial relic. 
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